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We present a theoretical analysis of several aspects of nonequilibirum cotunneling through a strong 
Coulomb-blockaded quantum dot (QD) subject to a finite magnetic field in the weak coupling limit. 
We carry this out by developing a generic quantum Heisenberg-Langevin equation approach leading 
to a set of Bloch dynamical equations which describe the nonequilibrium cotunneling in a convenient 
and compact way. These equations describe the time evolution of the spin variables of the QD 
explicitly in terms of the response and correlation functions of the free reservoir variables. This 
scheme not only provides analytical expressions for the relaxation and decoherence of the localized 
spin induced by cotunneling, but it also facilitates evaluations of the nonequilibrium magnetization, 
the charge current, and the spin current at arbitrary bias- voltage, magnetic field, and temperature. 
We find that all cotunneling events produce decoherence, but relaxation stems only from inelastic 
spin-flip cotunneling processes. Moreover, our specific calculations show that cotunneling processes 
involving electron transfer (both spin-flip and non-spin-flip) contribute to charge current, while spin- 
flip cotunneling processes are required to produce a net spin current in the asymmetric coupling 
case. We also point out that under the influence of a nonzero magnetic field, spin-flip cotunneling 
is an energy-consuming process requiring a sufficiently strong external bias-voltage for activation, 
explaining the behavior of differential conductance at low temperature: in particular, the splitting 
of the zero-bias anomaly in the charge current and a broad zero-magnitude "window" of differential 
conductance for the spin current near zero-bias-voltage. 
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I. INTRODUCTION 

Recent advances in probing and manipulating elec- 
tronic spin in semiconductor quantum dots and other 
nanostructures hold the promise of new applications re- 
lating to quantum computation and quantum informa- 
tion processing. A single electronic spin not only can be 
used as an elementary quantum memory unit, i.e., the 
qubit, due to its relatively long relaxation time in semi- 
conductors, but it is also expected to be useful as an ele- 
ment of calculation in the context of quantum computing 
algorithms or quantum information transport processing, 
which depend essentially on its temporal persistence of 
quantum interference^ 

This expectation provides strong motivation to develop 
a full understanding of the coherent evolution dynamics 
of a single spin in semiconductors. Actually, much ef- 
fort has been made on this matter from both theoreti- 
cal and experimental points of view. In particular, for 



a single quantum dot (QD), it has been predicted that 
measurements of the tunneling current between two leads 
via this QD may be an appropriate experimental tool to 
extract information about the orientation and dynamics 
of a single spin localized inside the QD. 2 Indeed, recent 
scanning tunneling experiments have observed modula- 
tion of the tunneling current through a single molecule 
with a spin subject to a constant magnetic field at the 
Larmor frequency, which is the characteristic dynamical 
(precession) frequency of a single spin under influence 
of magnetic field, with a corresponding peak in its noise 
power spectrum^* This feature has been examined the- 
oretically on the basis of two different weak tunneling 
models in the strong Coulomb blockade regime: (a) se- 
quential tunneling using a simple quantum rate equa- 
tion approach^ and (b) the two-channel Kondo Hamil- 
tonian using the nonequilibrium Green's function formal- 
ism jointly with the Major ana- fermion representation,*' 

However, from a fundamental quantum mechanical 
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point of view, any quantum measurement will inevitably 
introduce some disturbances into the measured system 
and consequently induce decoherence in the system vari- 
able conjugate to the one being measured. Therefore, the 
information concerning spin dynamics extracted from a 
tunneling measurement is expected to involve a reaction 
signature of the tunneling upon the coherent evolution 
of the single spin. It is therefore crucial to theoretically 
account for the tunneling-measurement-induced spin re- 
laxation and decoherence behaviors as functions of tem- 
perature and bias- voltage applied between the two leads; 
so as to provide a better understanding of the information 
obtained from measurement, and to give useful insight on 
how to raise measurement efficiency. 

The earlier papers cited above have concentrated on 
interpretation of the peak in the current noise spectrum, 
but a systematic investigation of the noncquilibrium re- 
laxation and decoherence effects, as far as we know, is 
still lacking. It is the main purpose of this paper to per- 
form these investigations. In treating the open quantum 
system at hand, we will employ the quantum Heisenberg- 
Langevin equation approach, to establish a set of quan- 
tum Bloch equations (i.e., equations of motion for the 
reduced density matrix) for a two-level system (a single 
spin in the QD) tunnel-coupled to two normal leads in 
a fully microscopic way, and then proceed to study the 
dynamics of a single spin qubit in an ambient magnetic 
field under nonequilibrium transport conditions. 

As mentioned above, two different tunneling mecha- 
nisms have been utilized to describe the quantum mea- 
surement process. If the chemical potentials of the two 
leads are nearly matched with the energy level of the 
sandwiched spin, the resonance condition is satisfied and 
the lowest-order tunneling, i.e., sequential tunneling, is 
observed in the transport process. However, it is quite 
likely that the chemical potentials are more probably far 
from the resonant point in actual experiments. In this 
case, the lowest-order tunneling of an electron into the 
QD is largely suppressed, but at very low temperature, a 
higher-order tunneling mechanism known as inelastic co- 
tunneling dominates the transport in the strong Coulomb 
blockade regime; in this mechanism an electron tunnels 
from the left lead to a virtual state in the dot, and then 
another electron tunnels from the dot to the right lead 
without changing the charge inside the QD. This is the 
tunneling mechanism that gives rise to the Kondo effect 
in QD tunneling. In fact, an exact mapping has been es- 
tablished between such cotunneling and the anisotropic 
Kondo problem by analyzing and comparing their respec- 
tive perturbation series for tunneling amplitudes^ Ac- 
cordingly, we will adopt the Kondo Hamiltonian in this 
paper to describe the inelastic cotunneling process and 
study its dissipation involved in coherent tunneling via 
the QD: we do so by developing a generic Langevin equa- 
tion approach in second-order perturbation theory with 
respect to the s-d exchange coupling constant, J, in the 
weak tunneling limit. 

Moreover, we will evaluate the nonequilibrium spin 



magnetization of a QD subject to a magnetic field in 
steady state and examine the behavior of charge flow 
(cotunneling current) through the QD within the same 
framework. Actually, an interesting calculation of the 
spin magnetization of a Kondo QD has already been 
carried out recently by means of the Majorana-fermion 
Green's function technique^ The striking result obtained 
in that paper is that their magnetization result differs 
from the thermal equilibrium formula even at zero or- 
der in the spin-leads exchange coupling, J. More im- 
pressively, theoretical analysis shows that if proper ac- 
count of this nonequilibrium magnetization is taken in a 
calculation of the current, the resulting differential con- 
ductance will demonstrate double peaks at bias voltages 
eV — ±gfiBB (the Zeeman energy), a signature of Kondo 
effect with a constant magnetic field B 7 even in the 
second-order perturbation calculations^ Of course, the 
log-signature peculiarity of the Kondo effect occurs es- 
sentially in the next orders of perturbation theory. Here, 
in the present paper, we ignore such higher-order terms; 
thus we confine our study to the ordinary cotunneling 
processes. A detailed analysis of the third-order pertur- 
bation contribution to the current has been established 
in Ref. ,9, computing an explicit logarithmic enhancement 
in current. 

Recent theoretical studies have shown that nonequilib- 
rium Kondo physics is fundamentally governed by weak- 
tunneling perturbation theory when the bias voltage is 
much larger than the Kondo temperature, Tk; this can be 
ascribed to current induced decoherence of the resonant 
spin- flip term in cotunneling processes, which can elim- 
inate the generic logarithmic divergence in conventional 
Kondo physics P I£L 1 ii2i This is the reason that the third- 
order contribution in Ref. provides a quantitatively rel- 
atively small modification to the second-order term in the 
nonequilibrium current formula in the weak tunneling 
limit (albeit qualitatively important). It is also shown 
in Ref. that the second-order calculation (cotunnel- 
ing) of the differential conductance exhibits cusps at bias 
voltages eV = ±j/jbB, a remnant of the Kondo effect. 
Therefore, we call this cotunneling behavior "Kondo- 
type" cotunneling. In this paper, we will systematically 
study this kind of nonequilibrium cotunneling through a 
single spin in a finite magnetic field using second-order 
perturbation theory and will specifically analyze the co- 
tunneling processes responsible for its special transport 
characteristics. 

In addition, we will also examine the behavior of the 
spin current and show that inelastic spin-flip cotunneling 
can produce a nonzero spin current for asymmetric cou- 
pling systems subject to a finite magnetic field and an 
activation bias-voltage. Unlike charge current, we find 
that the sign of spin current is independent of the direc- 
tion of the applied bias-voltage, but it does depend on 
the asymmetry of the coupling constants to the left and 
right leads and the direction of external magnetic field. 

The remaining parts of the paper are organized as fol- 
lows. In Sec. ^ we present the physical model used in 



3 



this paper: a single spin weakly tunnel- coupled to two 
normal leads. To focus on tunneling induced decoher- 
ence here, it is assumed that the single spin is free of 
any other dissipative heat bath except for the tunneling 
reservoirs. In Sec. lIIII we will then present the derivation 
of the quantum Langevin equations of motion for the sin- 
gle spin. In Sec. IIVI qualitative discussion and concrete 
calculations will be given for the resulting decoherence 
and relaxation rates as functions of magnetic field, bias- 
voltage, and temperature. Sec. \V\ focuses first on the 
derivation of closed-form expressions for charge current 
and spin current within the framework of the quantum 
Langevin equation approach developed here, and then 
addresses all possible cotunneling processes that occur in 
this system, and their respective contributions to the cur- 
rents. A numerical evaluation of differential conductance 
for the charge current and the spin current is provided in 
the last part of this section. Finally, our conclusions are 
summarized in Sec. IVII 



II. MODEL HAMILTONIAN 

We employ the two-lead Kondo Hamiltonian discussed 
above to model inelastic cotunneling through a single 
spin (or QD) in a magnetic field, B, in the weak-coupling 
regime: 

H = Ho + Hi, H\ = H Ie a + iftrans, (1) 

H xe fi= 22 ^VV [( c ^kT C i)k'T _ C vUl C nk'l)S 

+ C lk1 C v k' + C l ? ki C r ; k'T^' + ] ' 

iftrans = Jrl J^[( c flkT C Lk'T — c Rkl c Lk' l) & 
k,k' 

+ c Rk1 c Lk>l S ~ + c flkl c Lk'T 5+ ] + ( R ^ L )> 

where c^ kcr (c v ^ a ) creates (annihilates) an electron in lead 
i] (= L, R) with momentum k, spin a and bare energy 
£k- Jll, Jrr, and J LR = J RL = \J JllJrr are Kondo 
exchange coupling constants between the electrons and 
the localized spin-1/2, S = (S x , S y , S z ), S ± = S x ± iS y . 
Hq stands for the free Hamiltonian containing: (1) two 
noninteracting normal leads, individually in local equi- 
librium with temperature T (not to be confused with 
decay times to be introduced below) , respective chemical 
potentials ji^, and Fermi distribution functions defined 
as f v (e) = [1 + e^-^A's 7 "]- 1 ; and (2) Zeeman energy 
of the localized spin subject to magnetic field B (g and 
fj,B are the Lande factor and the Bohr magneton, respec- 
tively). It should be noted that we ignore the Zeeman 
effect in the lead electrons. The interaction part of the 
total Hamiltonian, Hj, also includes two terms: if re fl de- 
scribes the reflection processes, in which an electron from 
a given lead is scattered back into the same lead in both 



spin-conserving and spin-flip configurations; while if trans 
describes the transmission events, where an electron from 
one lead cotunnels into the other lead, also in both con- 
figurations. Except for the tunnel coupling, we ignore all 
other "environmental" decay interactions of the single 
spin. 

Here, we assume the leads to have a flat density of 
states, Prj, in the wide-band limit. We take the chemi- 
cal potentials, /i n , to vanish in equilibrium and use this 
choice as the reference of energy throughout the paper. 
In the nonequilibrium case, we assume the bias voltage is 
applied symmetrically, fiL = ~^r = eV/2. Throughout, 
we will use units with H = fee = e = 1. 

The conceptual structure of our model is predicated on 
the idea that the full system can be separated into two 
subsystems: one of which is the measured subsystem, the 
single spin, and the other consists of the two leads jointly 
comprising a "heat bath" or "reservoir" . The interaction 
between the two subsystems, ifi, must be weak in order 
that the separation of the two subsystems is physically 
meaningful. Accordingly, Hi generates dissipation for the 
dynamics of the "open" measured quantum subsystem, 
which is the principal focus for study in this paper. For 
notational brevity, we rewrite this term as a sum of three 
products of two variables: 



ifi = Q z F Q z + Q + F Q + + Q~F Q 



(2a) 



with 



Q* — Qrm' — ^w'( c ])kTS'k'T C rikl c ri'k' l) > 

(2b) 

Q + = QvV 1 = "W^kTS'k'l' (^c) 

r/,7/' ,k,k' 



as functions of reservoir variables, and the corresponding 
generalized forces, Fq, depend on the variables of the 
measured subsystem as 

Fq* — S z , Fq+ — S~ , Fq- = S + . (2e) 

Here, the terms Q ± Fq± describe spin-flip tunneling pro- 
cesses, in which the conduction electron spin changes its 
orientation in the process of tunneling, and the local- 
ized spin is also flipped. On the other hand, the term 
Q z Fqz is responsible for non-spin-flip (spin-conserving) 
tunneling, in which no spin exchange takes place. All of 
these tunneling processes are schematically elaborated in 
Fig. 3. 



III. QUANTUM LANGEVIN EQUATIONS 

In this section, we derive a generic quantum Langevin 
equation approach and establish a set of quantum 
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Bloch equations to describe the dynamics of a single 
spin modeled by Eq. Q. It is well known that the 
underlying quantum Langevin equation approach has 
been extensively developed and successfully employed 
in the contexts of quantum electrodynamics and quan- 
tum opticspiii^iii Albeit that the great advantage of this 
scheme is that it allows us to naturally incorporate the ef- 
fects of quantum noise introduced by the "environment" 
on the studied system variables^ we will take no ac- 
count of such fluctuation issues here. Considering that 
such noise has a very short correlation time (determined 
by the reservoir correlation time, r r ), it is reasonable to 
neglect it for the longer time scale (> r r ) of interest in 
the present paper. 

The Heisenberg equations of motion for the spin Pauli 
operators S z , and the lead operators are given by 



iS z = =: [S 2 ,i?i]_ : = : (Q~ S+ - Q+S' 

i5 ± = [S ± ,H]- =±A5 ± + : [^,#1]- : 
= ± AS ± ± : (2Q ± S Z - g z 5 ± ) : , 
ic^kT = [c v ui,H]- = e, ) kC,,kT+ : [c^kTj-Hi]- : 

— ^k^r/k^ 4~ ^ ^ (j?]j]C>gk'^ -)- e/^C^k't) 

k' 

+ S (i/jj^C^k'X + J \f] c f\V! x)] i 
iCjjkl = [ c r;kl,-ff]- = e7jkC,,kJ.+ : [Cnkl,Hi]- : 

— £r;k c r;kX "P ^ [ ^ (^r;r;Cr/k'X "I" JrjfjCfjk' J,) 



(3) 
(4) 

(5) 

(6) 



where we have e^k = £k ~ /■*«> ^ = -^(-R) if V — 

and A = g[isB. \A,B\± = AB ± BA are, respectively, 

the commutator and the anticommutator of operators 



A and B. The equations of motion for 



J ?7ka" 



are easily 



obtained by Hermitian conjugation of the equations for 
c 7)k<r- The colon-pair notation, : (• ■ ■ ) :, in these equa- 
tions denotes normal ordering of the operators, • • • , in- 
side the square brackets: all annihilation reservoir op- 
erators Cjjko- are placed to the right of all spin opera- 
tors, S z ^\ and the creation reservoir operators c^ kcr are 
placed to the left of all spin operators, if the operators in- 
volved have equal-time arguments. For instance, the last 
two lines in Eqs. JSJ and are already normal-ordered. 
This normal ordering employed here is an operator coun- 
terpart of determining a cumulant in terms of Feynman 
diagrams with the elimination of disconnected diagrams 
involving products of lower order Green's functions. The 
latter disconnected diagram terms involve the effects of 
weak coupling (to the bath) which oscillate rapidly at the 
high frequencies of microscopic dynamics, with attendant 
destructive interference. While such terms do contribute 
small quantum corrections ( "renormalization" , "radia- 
tive corrections" ) to the microscopic dynamics on a short 
time scale, they are negligible in the context of the much 
longer time scale implicitly under consideration in our 



formulation of a quantum Heisenberg-Langevin equation. 
A full explanation of the normal ordering scheme in the 
equations of motion is provided in Refs. ll3Hl7l to which 
we refer the reader. 

Formally integrating these Heisenberg equations of mo- 
tion from initial time to t we obtain the exact solutions 
for these operators as 



S z {t) =S z {0)-i / dt' : [S z (t'), ffi (*')]- : , (7a) 



S ± (t)=e* iAt S ± {Q)-i / dt'e^^*-^ 



c^a{t) =e^ kt Cj , kCT (0)-* / dt'e-^-^ 

Jo 

x : [^(0,^1 (*')]- : • 



(7b) 



(7c) 



In the absence of interaction, H\ — > 0, these solutions 
become: 



s* (t)=s z o(t% 

S±(t) = e ^ A{t - t ' ) St(t') 



C nk<r W 



-ie»,k(t— *')/-° 



r/k 



At')- 



(8b) 
(8c) 



A standard assumption in the derivation of a quantum 
Langevin equation is that the time scale of decay pro- 
cesses is much slower than that of free evolution, which is 
reasonable in the weak-tunneling approximation. In this 
context it is appropriate to substitute the time-dependent 
decoupled reservoir and spin operators of Eq. © into the 
formal solutions of Eq. JJJ . Obviously, the full solution 
for the reservoir operator comprises two contributions, 
one from free evolution and the other from reaction of 
the spin through the weak coupling, and we denote these 
with superscripts "o" and "i" , respectively: 



C?7ka" \t) 



(*) 



(<). 



with 



CT (t) = -i / dt' :[c nk<T (t),H I °(t')) 



(9a) 



(9b) 



where H° is composed of the operators in H\ which 
are replaced by their decoupled counterparts (interac- 
tion picture). In fact, this is just the operator formu- 
lation of linear response theory. It should also be noted 
that Eq. (|9a|l implies that the two subsystems, the quan- 
tum dot and the reservoirs, are completely isolated before 
<o = 0, and the perturbative interaction, Hi, is adiabat- 
ically switched on from the initial time t — to- Using 
Eq. ©, the reservoir variables, Q^\t), become (Ap- 
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(10a) 



(11a) 



5 ± (t) = St(t) + Sf(t) 



St(t)-i9(r) / dT:[S?{t),HfV)]-: 



(lib) 



pendix : 

Q^(t)=Q^(i) + Q^(i) 

= Q%, (t) - i0(r) y dr : [Q~, (t), fff(i')] 

Q±,(*)=0*) + Q^(i) 

= 0*) " T ^ : [Q^(i),W)]- :, 

J — oo 

(10b) 

with t — t — t' and 9(r) represents the Heaviside step- 
function. Similarly, the formal solutions for the spin op- 
erators are also divided into two parts: 

S'(t) = s*(t) + s?(t) 

= S z {t)-i9{r) f d T :[S' (t),HfW]-:, 



Substituting these approximate solutions of Eqs. i|10|) 
and (|llf> into the equations of motion for S z and S 
[Eqs. © and (0J] and taking average evaluations with 
respect to the reservoir electron ensemble (• • • ) e an d over 
the localized spin degrees of freedom (••■)«, one can de- 
rive the desired quantum Bloch equations up to second 
order in the Kondo coupling constant J. After some al- 
gebraic manipulations (details are provided in Appendix 
|A"|) . the quantum dynamic equations take the compact 
form: 



J 



(S z ) = -\e{r)f dr([Q-(t),Q+(t')} + ) e ([S+(t),F Q+ (t')U s 

z J -oo 

-\d{r) J^dT{[Q-(t),Qt(t')}-) e {[S + (t),F Q+ (t')] + ) s 
+\0{r) f dr([Q+(t),Q-(t')] + ) e ([S-(t),F Q - (t')].) s 

+ \0{t) J^dt'{[Qt(t),Q-(t')U e {[S-(t),F Q -(t')] + ) s , (12) 
(S ± ) = ^fiA{S ± ) =f 9(t) C dT{[Qt(t),Q*(t')] + ) e {[S*(t),F Q *(t')}-) s 

J — OO 

T 9(r) f dr([Qt(t),Qj(t')U e ([S z (t),F QT (t')] + ) s 

J — OO 

± l -9{r) £jr{[Ql{t),Ql{t')} + U[S ± (t),F Q 4t')}^ s 

± 1 -9{t) [ dT([Q^(t),Q^(0]-> e ([5 ± W,^(t')]+> s - (13) 



In these equations, we drop the superscript "o" in the 
spin operators occurring inside integrations, since they 
involve the dynamical spin variables after taking expec- 
tation values. However, it must be borne in mind that 
their time evolutions are governed by Eqs. (|8a|l and l|8bj) . 
Apart from free evolution, it is clear that the spin dynam- 
ics are modified by the spin-lead interaction in a way that 
is precisely relevant to the response function, R ab (t,t'), 
and correlation function, C ab (t, t'), (a, b = z,+, — ) of free 



reservoir variables, which are defined as: 

R a \t,t') = ±9(T){[Q a (t),Q b (t')]-) e , (14) 

C ab (t,t>) = ±9(r){[Q a (t),Q b (t')} + ) e . (15) 

These forms of quantum Langevin-type dynamic equa- 
tions, expressed explicitly in terms of the correlation and 
response functions of free reservoir variables, have also 
been proposed in Ref. by employing the quantum 
Furutsu-Novikov theorem. The present derivation seems 
more direct and its meaning is more transparent. 



6 



Considering the reservoirs to be in separate (local) 
equilibrium states except for differing chemical poten- 
tials (reflecting the bias- voltage driving the current) and 
noting that these free fermion reservoir operators, c° k(T , 

c rjk<r > obey Wick's theorem (without correlation between 
the leads), we can readily express the functions R ab (t, t') 
and C' ab (t,t') in terms of reservoir distribution functions. 
The calculational details are provided in Appendix [5] 
Here we cite some useful properties. Firstly, these re- 
sponse and correlation functions are functions only of 
the time difference r = t — t'. Secondly, these functions 
are related as 

R(t) = R+-(t) = R- + (t) = \r zz (t), (16) 

C(r) = C+-(r) = C- + (r) = \c"{r). (17) 

Therefore, it is convenient to introduce single Fourier 
time transforms for the two bath functions into frequency 
space: 



-4C(A) S Z + 2R(A). 



(22) 



/CO 
dre MT i?(r), 
-OO 

/oo 
dTe luJT C(r). 
-OG 



(18) 
(19) 



Thirdly, the spectral function C(u>) is an even func- 
tion of lo, while the imaginary part of the frequency- 
dependent retarded susceptibility R(u>) is an odd func- 
tion. In equilibrium, they are exactly related by the 
fluctuation-dissipation theorem. 

Employing the definitions of response and correla- 
tion functions and free evolution relation S' ± (t') = 
e ±lAr S ± (t), Eq. (E3 yields 

S z = -2S Z I dTe- lAT C(r) - f dre' lAT R(t) 

J — OO J — OG 

~2S Z I dre tAT C(T)+ ( dre lAT R(t). (20) 

J — oo J — oo 

(Hereafter, we suppress the brackets around the spin vari- 
ables since they are all c- numbers.) In a transport mea- 
surement experiment, a single spin decays to its exter- 
nal bias-voltage-driven steady state in a characteristic 
time, r c , of the system. If we assume that the single spin 
changes significantly only over a time scale f > r c , an 
appropriate Markov approximation may be generated by 
making the replacement 



dr 



dr. 



(21) 



From a physical point of view, this presumption is consis- 
tent with the normal ordering scheme performed in the 
operator equations of motion, in regard to elimination of 
the rapid oscillations in microscopic dynamics. In this 
case, the equation of motion for S z can be further sim- 
plified as 

S z = -2[C(A) + C(-A)]S Z + R(A)-R(-A) 



Analogously, the quantum Langevin equation for S be- 
comes 

-2S ± I dre ±tAT C(T) - 2S ± f d T C(r) 

J —OG J —OO 

= TiA5 ± -2[C(A) + C*(0)]5 ± . (23) 



IV. NONEQUILIBRIUM MAGNETIZATION, 
DECOHERENCE AND RELAXATION 

In this section, on the basis of the derived Bloch equa- 
tions, Eqs. (|22|l and 1)23(1 . we will carry out analyti- 
cal evaluations of the relaxation and decoherence rates, 
as well as the magnetization of the single spin under 
transport conditions, as functions of temperature, bias- 
voltage, and external magnetic field. 

It is well known that there are two decay mechanisms 
leading to standard Bloch equations which define two 
distinct relaxation time scales: (1) The longitudinal re- 
laxation time, Ti, is responsible for the spin magnetic 
moment relaxation, while (2) the transverse relaxation 
time, T2, is responsible for decoherence of the quantum 
superposition state composed of the two spin states a =t 
and J,. These time scales arc defined by the time evolu- 
tions of S z (t) and S^, respectively: 

-L = 4C(A) = 2vr ( j\ L p\ + J 2 RR pl) Tip (J± 
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V 



T 



LP 



V 



T 



,(24) 



±r = 2[C{A) + C(0)] = 4r + 2C(0) 
= ^ Fi +^{JIlpI + JIr.Pr)T 

2 fv 

+2tt J LR p L p R Tp> I — 



(25) 



In deriving these results, we employed Eqs. I(B19J| and 
(|B20p . It is noteworthy that the transverse spin re- 
laxation rate, I/T2, includes two contributions: the 
relaxation-induced dephasing, l/2Ti, and also pure de- 
coherence, 2C(0). 

From Eqs. (|T2l and ((T3l . we can easily deduce that 
the longitudinal relaxation time (Ti) stems completely 
from spin-flip cotunneling events, which is conceptually 
consistent with the physical definition of spin relaxation 
and implies its dependence on magnetic field. Further- 
more, spin-flip processes also contribute to decoherence 
with the partial rate, 1/2T%. In contrast, non-spin-flip 
processes do not induce spin relaxation but they do con- 
tribute to pure decoherence with the partial rate 2C(0), 
which is independent of magnetic field. This difference in 
the magnetic field dependences of the two rates may be 
understood in the following terms: the non-spin-flip pro- 
cess entails charge transport through the QD via a virtual 
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state but the QD eventually returns back to its original 
spin state without changing energy (which is why this 
process is referred to as elastic cotunneling in the litera- 
ture); whereas energy exchange does take place between 
the QD and leads in the spin-flip process in a finite mag- 
netic field, in which the spin of the QD is finally flipped 
and thus the QD is inelastically excited or decays ac- 
companied by excess energy, the Zeeman energy, A. Of 
course, in the absence of an external magnetic field, spin- 
flip cotunneling also becomes elastic. In this case, the two 
relaxation times are equal (i.e., it is inelastic cotunneling 
that makes them differ), 



4C(0) = At: (j 2 LL p 2 L + J 2 RR p 2 R ) T 

'V 



+4nJ LR p L PRT(p 



T 



(26) 



On the other hand, in the limit of zero bias-voltage 
(equilibrium condition), the relaxation rates become 



f 

TP 



= 47T 



JllPl + JrrPr 



JlrPlPr ) T(p [ — 



(27) 



rpcq 
1 2 



2T cq 



4 JM+« + J 2 LR p LPR \ T. 



(28) 

It is clear that the contribution of non-spin-flip cotun- 
neling to the pure decoherence rate, the second term in 
Eq. i|28|) . is proportional to temperature, while spin-flip 



cotunneling leads to a somewhat complicated tempera- 
ture dependence, Acoth(A/2T). This difference can also 
be ascribed to energy exchange in the dissipation process. 
It is worth noting that without transport (V — > 0), dissi- 
pation (relaxation and decoherence) is due solely to quan- 
tum thermal fluctuations (thermal noise). Furthermore, 
if the external magnetic field is quenched, the thermal 
fluctuations are purely non-energy-consuming: 



Tth 



_L = ±_ = 8lT ( JLpI + JrrpI + ji rPlPr ) Tj 



2? 



indicating that the dissipation is totally determined by 
thermodynamics (temperature) of the reservoirs. 

For illustrative purposes, we exhibit in Fig. 1 the 
dependences of the relaxation rate and the dephasing 
rate on magnetic field (a,b) and on bias-voltage (c,d) 
for given temperatures. At relatively low temperatures, 
T/V = 0.01 (or T/A = 0.01), these rates show linear in- 
crease with respect to bias-voltage V (magnetic field A) 
with the rates of increase depending on the relative mag- 
nitudes of V and A. Interestingly, dephasing, 1/T 2 , is 
independent of V for V < A, as shown in Fig. I (c). This 
comes about because the hyperbolic cotangent functions 
behave as <p((A + V)/T) + tp((A - V)/T) -> |A + V\ + 
\A — V\ in the limit T — > 0. As expected, rising tem- 
perature smears out the low-temperature structures in 
these rates [Figs. I (b) and (c)]. Finally, the temperature 
dependences of the two rates are summarized in Fig. 2. 

The magnetization of the QD, defined as M — S z , is 
readily obtained using the steady solution of Eq. lf^2l as 



M(A,V)=S z oc = 



g(A) 
2C(A) 



J LlPl + JrrPr t2 \ A 

2 RR R + JIrPlPr) 



{JllPl + j2 rrPr) V 



T 



JlrPlPr 



V 



T 



■ip 



V 



T 



(30) 



which is identical to previous theoretical result&S In ab- = i((Q^ R — QrI)S z ~ (Qll + Qrl)^ + 

sence of bias- voltage, V = 0, it reduces to the equilibrium +{Qll + Qlr)S~)^ (31) 

expression M(A,0) = |tanh(A/2T). 



V. NONLINEAR TUNNELING CURRENT 

The calculation of steady state tunneling charge cur- 
rent, I c , measuring the charge flow from left lead to right 
lead, is based on the equation of motion for the charge 

density N L = Y.a N Lv = Ek c Lk/Lk«r in the left lead > 
which is a sum of both spin-up and spin-down electrons 
flows, I c = Tlt + Ilu 

/l T = -(N^)=i([N Lh H}) 



Ili = -(N Ll )=i([N H ,H]) 

= - Q z lk)s z + (Ql L + q- lr )s + 
-(Q + ll + Q + r L )s-), (32) 

with the definition Q*°? = J w £ k k , c\ Va c , kV . Using 
the same procedure described in the preceding section 
and Appendix ^ and employing the various response 
and correlation functions of the free reservoir variables 
determined in Appendix El we have derived analytic ex- 
pressions for the spin-resolved current, 1^. For example, 
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FIG. 1: The cotunneling-induced spin relaxation rate (T-^ 1 ) 
and decoherence rate (T^ - ) as functions of magnetic field A 
(a,b) and of bias- voltage V (c,d) for given temperatures indi- 
cated in these figures. The parameters we use in calculation 
are: JllPl = JrrPr = Jlr^PlPr = 0.02. Here, we use the 
decay rate of purely thermal fluctuations, r t h, Eq. (J2HJ, as the 
unit of the two rates. 



the spin-up current, takes the form 



dr 



(It 



LR,RL\'J ^RL.LRVJ 



R-+ 
n LL,LL 



(r) + RrIlr^) 



-e 
-2S 



(At 



} 



(At 



c 



LL,LL 



{e' lAT [c- L t LL {r)+C R l LR {T) 
(t)+C+ r . rl (t)]}. 



(33) 



Using Eqs. (|B4|) . (|B5|I . (JB9J (JBTTJ , and then making the 
replacement J^^dr =>■ dr, we finally arrived at an 
explicit result for 1^ as a function of temperature and 
bias-votage [after performing the e-integrals of Eqs. 1|B16|) 
and JEHJl], 



3tt 



'LI 



2 

2ttT 



JlrPlPrY + ttA (Jii/s! + J 2 lr PlPr) 



JllPlV I y 



JlrPlPrV 



A + V 
T 



S z . 
(34) 



Similarly, the current from spin-down electrons, is 
given by 



37T 

hi = y J lrPlPR v ~ 7:A ( j llPl + JlrPlPr) 

t2 2 /^ A \ ,2 ( A~V 

JllPlV j lrPlPr<P ( T 



2nT 



S z . 
(35) 



The total charge current, I c , is thus 



1° = nJl R p L pR<3V + 2S z T 



A-V 
T 



A + V 
T 



(36) 



It is worth noting that the cotunneling current, Eq. 13011 . 
is just proportional to second-order in the exchange cou- 
pling constant, Jlr, and higher-order contributions are 
all neglected. This is because we employ the approxima- 
tion formula, Eq. (|A2|) . to derive non-Markovian quan- 
tum dynamic equations and the current, leading to the 
absence of the characteristic logarithmic divergence term 
in current. Thus, the present approach can not be ap- 
plied to describe strong Kondo correlations, but can be 
used to study the ordinary cotunneling process in the 
weak tunnel-coupling limit. 

Due to the fact that spin-up electrons are coupled with 
spin-down electrons via the spin-flip processes in this 
model, there is an imbalance between the spin-up current 
and spin-down current, i.e., there is a net spin current, 
I s , with respect to the left lead, defined as 



I s = Il\ - hi = 2ttA{JI lP 2 l + JIrPlPr) 
~2ttT I 2 Jl L p 2 L ip ( — j + JIrPlPr V 



'A-V 
T 
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-<P 



A + V 
T 



S z 



(37) 



To better understand these formulae and the physical 
perspective involved, we elaborate the physical picture 
of cotunneling processes through a QD in a finite mag- 
netic field. When the electronic levels of the QD, e<j CT , 
are far below the chemical potentials of the two leads, 
i.e. tda *C Pl(r)i the first-order tunneling process, se- 
quential tunneling, vanishes. However, higher-order tun- 
neling processes, cotunneling, are active and dominate 
the quantum transport. In the strong Coulomb blockade 
regime, the QD is always singly occupied by an electron 
because of the deep electronic energy, e^, and the ex- 
tremely strong charging energy, U — > oo, involved when 
an additional, excess electron attempts to enter the QD, 
i.e. €da + U P-l(R)- This steady occupation means 
that no charge-fluctuation takes place, provided that the 
applied bias-voltage is not strong enough to force the 
chemical potential in one of leads below the QD level, 
tda-, so as to drive the transport into the sequential tun- 
neling regime. Therefore, in nonequilibrium conditions, 
a cotunneling event consists of two single-particle tun- 
neling processes, (Tj and @, which can take place in se- 
quence as follows: event @, an electron inside the QD 
with spin-er will at first tunnel out to a lead L or R, 
inducing a virtual empty state in the dot, but this is im- 
mediately followed by a second single-particle tunneling 
event (2) in which an electron in one of the leads is in- 
jected into the QD with the same spin-er (spin-conserving 
elastic cotunneling) or with the opposite spin (spin-flip 
inelastic cotunneling). The two tunneling events occur 
via a virtual empty-dot-state in a very short time inter- 
val to insure coherence. Importantly, spin-flip cotunnel- 
ing provides a mechanism for the spin orientation of the 
QD to be changed, which is fully quantum phenomenol- 
ogy. Moreover, it stimulates intrinsic spin-fluctuation, 
which is a fundamental concept in the context of Kondo 
physics. Obviously, there is a total of 16 different cotun- 
neling events allowed in this strong Coulomb blockade 
system, which are schematically shown in Fig. 3. We can 
classify them in four different categories/types. 

Type-I cotunneling involves only one lead and equal 
spin orientations in two successive single-particle tunnel- 
ing events. Figures 3(a)-(d) depict such trivial cotunnel- 
ing events, in which an electron with spin-cr exits the QD 
to the lead r\ and subsequently an electron with the same 
spin in the same lead (probably not the same electron) 
transfers to the QD. Obviously, these events make no con- 
tribution to the current. In contrast, type-II cotunneling 
experiences a spin-flip process as shown in Fig. 3(e)-(h), 
in which the final spin state in the QD is opposite to its 
initial state. Of course, only events (e) and (f) relate 
to the currents in the left lead, and involve the terms 
TQl L S + and ±Qt L S~ in Eqs. (SJ and However, 
the spin-up and spin-down electrons flow in opposite di- 
rections and contribute to the currents with equal magni- 
tudes. As a result, no charge current occurs in this type 
cotunneling but spin current does emerge. Furthermore, 



because only one lead is involved, the contributions of 
events (e) and (f) are naturally independent of the bias- 
voltage and are only dependent on the magnetic field, 
involving the terms ±TxAJ\ L p\ qp ^TJl L p 2 L tp(A/T)S z 
in spin-up and spin-down currents with ± — * +, — , re- 
spectively. Figures 3(i)-(l) represent all type-Ill cotun- 
neling processes. This kind of cotunneling describes an 
equivalent spin-conservative tunneling process, in which 
an electron is transferred from one lead to another lead 
via the QD without spin exchange. The corresponding 



ztT(U) 

LR(RL) 



S z . More- 



terms in Eqs. (|3*T)l and are ±Q 

over, we observe that (1) the spin-up events (i)&(j) and 
the spin-down events (k)&(l) yield currents having not 
only the same directions, but they also have equal mag- 
nitudes; and (2) these spin-up(down) events involve only 
the difference of chemical potentials, leading to a contri- 
bution proportional to the bias-voltage, ^TTj 2 R pLpfjV. 
This observation reveals that type-Ill cotunneling ex- 
cludes the possibility of spin current, but it does provide 
a linearly bias- voltage-dependent term in the charge cur- 
rent, Eq. 1)36(1 . Finally, type-IV cotunneling comprises 
the four electron-transferring tunneling events accompa- 
nied by a spin-flip process as exhibited in Figs. 3(m)-(p). 
They produce the terms —Q R ' L S + , +Qt R S~ in Eq. IpTTjl 
and QZrS + , ~Qrl^~ m Ec l- (E2J - Differing from type- 
Ill cotunneling, we find that in type-IV cotunneling the 
spin-up(down) events (m)&(n) [or (o)&(p)] involve both 
the voltage change and spin-flip, and the corresponding 
contributions to current are dependent on both V and A: 
±[A - 2S z Ttp(A ± V/T)]nJl R p L p R . Type-IV cotunnel- 
ing produces both spin and charge currents. In sum, the 
mechanism for creating spin current stems solely from 
inelastic spin- flip cotunneling processes (type-II and IV), 
while the electron-transferring elastic and inelastic co- 
tunneling processes (type-Ill and IV) are responsible for 
producing charge current. 

Substituting the steady-state solution, Eq. (|30|l . into 
the charge current, Eq. I|36|) . and the spin current, 
Eq. (|37|) . we readily find that: (1) both the charge cur- 
rent and the spin current are zero if V = 0; (2) the 
resulting spin current is nonzero in nonequilibrium con- 
ditions, V 7^ 0, if and only if two conditions are satisfied: 
Jll 7^ Jrr, i-e. the asymmetrical Kondo coupling case, 
and there is a nonvanishing magnetic field, A ^ 0; (3) 
the spin current is an even function of the applied bias- 
voltage, V , indicating that the sign of spin current is 
not related to the direction of the bias-voltage; whereas 
charge current, Eq. 1(36(1 . is an odd function of the bias- 
voltage and will change its sign when bias-voltage is ap- 
plied in the opposite direction; (4) the magnetic-field- 
related spin current changes its sign when the direction 
of the applied magnetic field is reversed (odd function) , 
while the charge current is an even function of A, because 
it measures the results of the total charge flow irrespec- 
tive of the spin orientation. The sign property of the spin 
current was also pointed out in previous study^S 

As an illustration, we plot the bias-voltage dependent 
differential conductance, dI c /dV, in Fig. 4. The differen- 
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FIG. 3: Schematic description of all 16 cotunneling processes (a)-(p) through a QD subject to Zeeman splitting energy A 
between the spin-up and -down electronic states, in the strong Coulomb blockade regime. A finite bias-voltage V is applied 
between the two leads, L and R. An open circle inside the QD stands for the initial occupied electron state before tunneling 
events, while the solid circle denotes the final state occupied by an electron after the cotunneling processes. The single-particle 
tunneling event (T) takes place first and then is followed by the tunneling event (2), which together comprise the entire cotunneling 
process. The reservoir variable below each of the figures denotes the corresponding physical process occurring in the reservoir. 
The arrow beneath the left lead in each of the figures denotes the flow direction of the spin-up and(or) spin-down electron with 
respect to the left lead. 



tial conductance shows a characteristic jump at V — ±A, 
which is the signature of the Kondo effect in the pres- 
ence of an external magnetic field. Mathematically, this 
feature comes from the hyperbolic cotangent function in 
the current formula, Eq. (|36fl . From a physical point 



of view, this splitting can be qualitatively understood 
from the following consideration: a small bias-voltage, 
|V| < A, can not provide enough energy to spur the 
spin-flip cotunneling process that is an energy-consuming 
event in the case of nonzero magnetic field; However, 
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when \V\ > A, the spin-flip cotunneling process is ener- 
getically activated, thus an additional channel is opened 
for electron transport. Moreover, the effect of tempera- 




FIG. 4: The calculated differential conductance dI c /dV vs. 
bias-voltage V/A for several temperatures at nonzero mag- 
netic field in units of Go = ^ttJlrPlPr (linear conductance 
at zero magnetic field). Other parameters are the same as in 
Fig. 1. 
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FIG. 5: The calculated spin current, I s , and its differential 
conductance, dl a /dV, as functions of bias-voltage V/A, at 
nonzero magnetic field, (a) exhibits results for several tem- 
peratures and Jrr/Jll = 4.0, Jll = 0.02. (b) plots the 
results for Jrr/Jll = 5.0 (Jll = 0.02) as solid lines, and for 
Jll/Jrr = 5.0 (Jrr — 0.02) as dashed lines. 

We also exhibit the resulting spin current, I s , and its 
differential conductance, defined as dl s /dV, as functions 
of bias- voltage V/A in Fig. 5. We observe that at low 
temperatures, T/A = 0.01 and 0.05 in Fig. 5(a), the 
calculated spin currents are nearly zero in the small bias- 
voltage region, \ V\ < A, notwithstanding Jll/Jrr = 4.0 



and A 7^ 0. Analogous to the peak-splitting of the differ- 
ential conductance shown in Fig. 4, the low-temperature 
vanishing of spin current is also due to the fact that spin- 
flip scattering is energetically inaccessible in the case of 
small bias-voltage. This vanishing produces a "window" 
of zero differential conductance for spin current. Never- 
theless, at higher temperatures, thermal fluctuation pro- 
vides an additional possibility to flip spin in the tunneling 
processes, leading to a slow increase of the spin current 
and the gradual disappearance of the zero "window" in 
dI s /dV. In Fig. 5(b), we show that the sign of the spin 
current is also determined by the relative magnitudes of 
Jll and Jrr. 



VI. CONCLUSIONS 

In this paper, we have systematically examined 
nonequilibrium inelastic cotunneling through a single 
spin (QD) subject to a finite magnetic field in the strong 
Coulomb blockade regime, in the weak tunnel coupling 
limit. For this purpose, we introduced the Kondo Hamil- 
tonian to model cotunneling in the QD and employed 
a generic Heisenberg-Langevin equation approach to es- 
tablish a set of quantum Bloch-type dynamical equations 
describing inelastic cotunneling phenomenology. 

In our formulation, the operators of the localized spin 
and the reservoirs were first determined formally by inte- 
gration of their Heisenberg equations of motion, exactly 
to all orders in the tunnel coupling constants. Next, 
under the assumption that the time scale of the decay 
processes is much slower than that of free evolution, we 
expressed the time-dependent operators involved in the 
integrands of these equations of motion approximately in 
terms of their free evolution. Thirdly, these equations of 
motion were expanded in powers of the tunnel- coupling 
constants to second order; this approximation is physi- 
cally valid in the weak tunnel-coupling limit. On the ba- 
sis of these consideration, jointly with normal ordering, 
we developed the Bloch-type equations expressed explic- 
itly and compactly in terms of the response and the cor- 
relation functions of the free reservoir variables, which 
facilitated our theoretical examination of relaxation and 
decoherence in the localized spin induced by the "envi- 
ronment" . 

In the problem at hand, dissipation of the QD spin 
stems from tunnel-coupling of the QD to two leads by 
cotunneling mechanisms. Based on our derived Bloch 
equations, we obtained explicit analytical expressions for 
the corresponding relaxation and decoherence rates at 
arbitrary bias-voltage and temperature. We found that 
relaxation results exclusively from spin-flip cotunneling 
processes alone, whereas both spin-flip and non-spin-flip 
cotunneling events contribute to decoherence. In this 
analysis, we carried out systematic examinations of the 
relaxation rate and the decoherence rate as functions of 
bias-voltage, external magnetic field, and temperature. 
Our formulation also facilitated the derivation of an an- 
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alytic expression for the nonequilibrium magnetization 
that is found to match that of earlier theories. 

Employing this approach, we also derived closed-form 
expressions for the spin-resolved currents, which facil- 
itated our calculation of both the charge current and 
the spin current. Furthermore, we classified and exam- 
ined all possible cotunncling processes occurring in the 
strong Coulomb interaction QD (16 events), and cate- 
gorized them in four distinct types. In this, we found 
that: (1) type-I cotunneling make no contribution to 
current; (2) spin-flip processes, types-II and IV cotunnel- 
ing, drive the spin current; (3) the electron-transferring 
processes, types-Ill (non-spin-flip) and IV (spin-flip) co- 
tunncling, produce charge current; and (4) we also de- 
termined formulae for their respective contributions to 
the charge and spin currents. Our numerical calcula- 
tions exhibit splitting of the zero-bias-voltage peak in 
the differential conductance for charge current in a finite 
magnetic field, which is a typical signature of the Kondo 
effect, and a wide "window" of zero differential conduc- 
tance for spin current about zero-bias- voltage. With in- 
sight gained from our specific analyses, we can ascribe 
these low-temperature transport characteristics to the 
fact that inelastic spin-flip cotunneling is energetically 
active only for sufficiently strong applied bias-voltage, 
V > A. We have also shown that spin current, unlike 
charge current, is an even function of the applied bias- 
voltage, and its direction depends on the orientation of 
the ambient magnetic field and asymmetry of the Kondo 
coupling constants to the left and the right leads. 
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APPENDIX A: DERIVATION OF EQS. itTOal) . O, 
AND lfl3|l 

In this appendix, we first prove Eq. (|10fl . Consider 
Q z m , for example. Substituting Eq. JJJJ into the definition 
of Eq. (|2b|) , we have 

Q z m -{t) = JW£[4t + 

k,k' 

x[^ T (t) + ^tW]--WE(H) 

k,k' 

= <C(t)-iJWE / dt ' 

k ,k' J ° 

- ■ W)]- : c£ k , T (t)} +iJW$>-i) 

k,k' 



= Q%(t)-iJ vv , y £9(T) dt' 

k,k< J ° 

x : [c&tW^vtW,^*')]- : 

k,k' 

= Q-,(i) - *0(r) f dr : [Q%,(t),Hf{H)\_ : . 

J — OO 

(Al) 

In the second stage of Eq. (|A1|) , we neglect terms of the 
form, c^ kcr cj ) , k , cr , , since they are of second-order in the 
coupling constant, 0(J 2 ), yielding a third-order contri- 
bution to Q* , with respect to J. 

To derive Eqs. (|12f> and (|13l) . we consider the nor- 
mally ordered product of the reservoir and spin opera- 
tors, : Q a S b :, 

: Q a (t)S b (t) :=: [Q«(t) + Q?(tM(t) + : 
= :Q a (t)S b (t) :-i6(T) [ dr {: Q a a (t) 

J —OG 

x : [s b Mzm]- ■■ + ■■ mim{t')]- ■■ ■.} . 

(A2) 

Once again, we neglect the term QfS b as it is propor- 
tional to 0( J 3 ). The first term in Eq. I)A2|1 involves only 
the free reservoir variables and the decoupled single spin. 
The other interaction terms (we designate the operator 
expressions in the integrand as X) arise from tunneling 
reaction, upon which we focus in the following derivation. 
Using the compact definition of the interaction Hamilto- 
nian, Eq. (J2EJ) , Hi = J2ce{z.+.-} Q CF Q C > we nave 

^ = J2{--Q a o(t)Q c o(t')[s b (t),F^(f)}_ 

+ [Q a (t),Q c (t')]^ c (t')S b (t):} 

= Ef \[Qi{t)mt')] + [s h {t),F° Qe {t')]^ 

c 

+ \[Ql(t)mt')]-[S h {t),F^{t')\- 
+ [Q a (t),Q c (t')UF^(t')S b (t):} 

= ^El : [Qo(t),Q c (t')} + [s b (t),F^(t')}^ 

MQ a o(t),Q L (t')}-[S b (t),F^(t')} + :} . (A3) 

Therefore, the full normal-ordered operator product 
(: Q a S b :) is written as the sum of a zero-order term 
and a term of second-order in the coupling constant J, 
having the compact form: 

: Q a (t)S b (t) :=: Q a (t)S b (t) : -i [ dr^i- C ac (t,t') 

J-oo c 

x[S b (t),F^(t')]- + R ac (t,t')[S b (t),F£ c (t')] + :}, (A4) 
with the definitions 

R ac (t,t<) = ±9(r)[Q a (t),Q c (t')U (A5) 
C ac {t,t') = \e{r)[Q a {t),Ql{t')] + . (A6) 
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The reservoir equilibrium ensemble averages of R ac (t,t') 
and C ac (t, t') are just the response function R ac (t, t') and 
the correlation function C ac (t, t') defined in Eqs. 114|l and 
(fTo*)) . respectively. 

In the spin operator equations of motion, these zero- 
order terms contribute quantum fluctuations associated 
with the reservoir fields, as well as quantum effects per- 
taining to the intrinsic character of the reservoirs (for 
example, superconducting or ferromagnetic leads). In 
any event, we take the ensemble average of each equa- 
tion of motion separately in regard to the electron en- 
sembles of the reservoirs and in regard to the quantum 
spin states. Thus, the normally ordered operator prod- 
ucts factorize in the averaging procedure. Considering 
that we take no account of quantum fluctuations in the 
present paper and only normal leads are connected to 
the single spin, the zero-order terms make no contribu- 
tion to the quantum Bloch equations. Moreover, only 
(Qt(t)Q*(t')) e and (Q*{t)Q z (t')) e are nonzero for nor- 
mal leads (see Appendix [BJ) . Combining all the above 
results, we obtain Eqs. i|12fl and (|13|1 . 



APPENDIX B: RESPONSE AND CORRELATION 
FUNCTIONS OF THE RESERVOIRS 



To obtain explicit expressions for the Bloch equations 
and the current, we need to determine the various cor- 
relation and response functions C(R)^„ 2 „ 4 (t, t 1 ) of the 
free reservoir variables, which are defined as 



(Bl) 

In the following, we drop all super(sub)scripts, "o", 
bearing in mind that all operators are free reservoir op- 
erators. In our calculations, we assume that: (i) the leads 
have a flat density of states p v for both spin orientations, 
so we can make the replacement 



£(•••)—/>., / &(•••); 



(B2) 



(ii) the normal leads are in the respective bias-voltage 
driven local equilibrium states described by 



m = 



l _|_ e (e-/*„)/T 



(B3) 



with temperature T and chemical potential fj, v ; and (iii) 
the time evolution of free reservoir operators is governed 
by Eq. (|5cl) . According to Wick's theorem and proper- 
ties (i) and (ii), it is easy to see that only the functions 
C(-ft)^',77'7)! C(R) n y , v , and C(R) m ,^, n 
We calculate them individually. 



are nonzero. 



i. c(Ry*(t,t') = j:^ / c(Ry ri *,, ri , rl (t,t'y. 
c(R)T£ RL (t,t') = \e{r){[Q^ R {t)M^i{t')]±)e 

= \0{r)J 2 LR J2 ([ c ikTW^ T W' c kT^)c iqn (i')]±> 

k,k ; ,q,q' 

= \0{r)J 2 LR e^«'-^ 

k,k' ,q,q' 

x [( c k T W c Lq< T (*))« ( c wt W c kt (*))« 

± (4qT M C flk'T CO) e ( c Lq' T (*)4kT W) e 

= \0{t)JIrPl P r J dede'e^> 
x{f L (e) [l-f R (e')]±f R (e')[l-f L (e)}}. 



Exchanging the roles of R and L, L 

)zz\\ , 
RL,LR\ 



C(Ry R z P LR (t,t') yields 



(B4) 
R, 



C{R)%g LR (tJ) = \e{r)Jl RPLPR Jdede'e-^-<> 
x {/ fl (e') [1 - Me)] ± f L {e) [1 - f R (e')}} . (B5) 

As we take the leads to be normal met- 
als/semiconductors, we have 



C{RT L % L {t,t') = \e(r)([Ql^(t),Q R ^(t')] ± ) t 
= C(R) z L z ^ RL (t,t'), 

)zz\\ , 
R L,L R \ 



c{R) R t l LR {t,t') = C(R) R ^ LR (t,t>), 



R L,L R \ 

C(R& RL (t,t') 



>L R , R L 
L R , R L 

Furthermore, 



C(R& RL (t,t') 



0. 



T){mmM R [ t{t')]±)e = o 



C(R)lU L (t,t>) = -d(r)([Ql L (t),Ql L (t')] ± ) e 
= <?(r) Jl L pl J dede'e^> 

x {f L (e) [1 - f L {e')] ± f L (e>) [1 - Me)]} , (B6) 



CWWM') = f(r)([Q z RR (t),Q RR (t')] ± ) e 
= 0(r)J 2 RR p 2 R J dede'e^'^ 

x {f R (e) [1 - Me')] ± Me') [1 - Me)]} ■ (B7) 

Finally, C(R) ZZ (t,t') are functions only of the time dif- 
ference r and take the form: 

C{Ry*{T)=e{r)Y,J 2 m Pl J dede'e^-^ 
x{/n(e) [l-/,(e')]±/,(e')[l-/,(e)]} 
+0{r)Jl R PLP R I dede' \e^> ± e -*(«-«')r 
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x {f L (e) [1 - / fl (e')] ± MO [1 - Me)]} . (B8) 
2. C(ii)+-/-+(t,f) = C{R)^{t,1f): 

C(R)1r&(t) = \e{r){[Q + L '-{t),Q-J L + {t')] ± ) e 
= \0{t)JI R PlPr J dede'e^'^ 
x {Me) [1 - MO] ± MO [1 - Me)]} > (B9) 



x {/ fl (e) [1 - MO] ± MO [1 - M*)]} . (BIO) 



C(i?)+^ + (r) - ^(r)([Qir(t),Qir(0]±)e 

RR,RR Z RR RR 

= \0{r)J 2 LL/RR pl /R f dede'e^'^ 

x {Mfl(e) [1 - /z/b(0] ± Mb(0 [1 - /z/h(<0] } • 

(Bll) 

Moreover, we can easily obtain 



)+/- 



7+^ 



= ^E -OS / & [A( e ) - ^ + w )] 

+ ^ 7r JlrPlPr J de[f R (e - u) - f R (e + u)] . 

(B14) 



It is readily seen that 



C(-w) = C(cj), fl(-w) = -i?(w). 



Using the formulae 



(B15) 



and 



de[f v (e + w) - f v >(e)] = -(/v - /u„ + w), (B16) 



+ <*)[!- Me)] = (BIT) 



we can perform the e-integrals in the large bandwidth 
limit, with the results 



C(R) + -(t) = C(R)-+(t) = -C(R)"(t). (B12) 

Therefore, in w-Fourier space, the spectral function 
C(ui) defined in Eq. (fT7)l is given by 

C(uj) = lf dTe lUT C zz {T) 

^ J — oo 

= \* E J lnP% I de iMe) [1 - + w)] 

+ /,(€ + W) [l"/ n (c)]} 

+ 2~* J lrPlPr J tfeE{/Tj( e )[ 1 -/*?( e + w )] 



+/„(e)[l-/ fl (e-w)]}. 



(B13) 



Also, the imaginary part of the retarded susceptibility 
R(uj) defined in Eq. JTSJ) is 



R(u) = ~ [ dre luJT R zz (T) 

2 J-oo 



fl( W ) 

CM 



;r| ' ^p|yvl + JLp£Pfl)u; > (BLS) 



| (JIlPl + JrrPr) Tip 



-JlrPlPrT 



■•P 



uj + V 
T 



<P 



uj-V 



T 



(B19) 



where we have defined 



tp(x) = aicoth . (B20) 



It should be noted that when the leads are in thermody- 
namic equilibrium, [i^ = PR, the spectral function C(lo) 
and the function R(u>) obey the fluctuation-dissipation 
theorem C{w) = R(u) coth(w/2T). 
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